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ABSTRACT 

We have included opacity tables in our stellar evolution code that enable us to ac- 
curately model the structure of stars composed of mixtures with carbon and oxygen 
independently enhanced relative to solar. We present tests to demonstrate the effects 
of the new tables. Two of these are practical examples, the effect on the evolution of 
a thermally pulsing asymptotic giant branch star and a Wolf-Rayet Star. The changes 
are small but perceptible. 
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1 INTRODUCTION. 

The main source of a star's energy is the nuclear fusion re- 
actions occurring either in its core or in thin burning shells 
around the core. This energy is transported from the produc- 
tion site to the surface by radiative transfer or convection. 
The first transports energy in the form of photons, while the 
second is a cyclic macroscopic mass motion that carries the 
energy in bulk. 

In regions where convection is stabilised, radiative 
transfer leads to the equation of stellar structure, 

dT _ 3krp L r . 
~dr ~ ~ 16aT s A-xr 2 ' ( ' 

where L r is the luminosity at radius r, p the density, a 
the Stefan-Boltzmann constant, T the temperature and Tir 
is the Rosseland mean opacity. 

Opacity k is a measure of the degree to which matter 
absorbs photons. There are four main sources of opacity in 
stars 

• Bound-bound transitions are the transitions of an elec- 
tron in an atom, ion or molecule between energy levels which 
are accompanied by either the absorption or emission of 
a photon. Only a photon with the correct wavelength can 
cause a given transition, the process is wavelength depen- 
dent. 

• Bound-free transitions, or photoionisation, occur when 
an incoming photon has enough energy to ionise an atom or 
ion and free an electron. The reverse process is the capture 
of an electron by an atom or ion. It will not occur until 
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photons above a threshold energy are available and falls off 
as i/~ 3 where v is the frequency. 

• Free-free transitions are scattering processes which oc- 
cur when an electron and photon interact near an atom or 
ion. The process is also known as bremsstrahlung. Again it 
is proportional to v~ A . 

• Electron scattering is wavelength independent at low 
temperatures where it is Thomson scattering. The electron, 
with its low cross-section, is a small target and so only dom- 
inates at high temperatures when most atoms are ionised. 
At very high temperatures relativistic effects are important 
and Compton scattering dominates. 

From this list it is possible to see that calculating the 
opacity in a stellar model is a difficult process and depends 
on the composition, temperature and density of the material. 
Deep in the star local thermodynamic equilibrium (LTE) is 
achieved and an average for radiative transfer ov er all wave- 
lengt hs requires the Rosseland mean opacity kr ( R.o sselandl 
1924) expressed as 

r i dB 

j_ _ I ku dT v 

Jo dT V 

where k„ is the opacity at frequency u, T is the temper- 
ature and B v is the flux per unit area into unit solid angle 
per unit frequency in LTE. 

The most c omprehensive opacities availab le today are 
from the OPAL Jlglesias fc Rogersll996t) or OP dCunto et all 
1993) groups who have made detailed models of the above 
processes. They provide tables of the Rosseland mean opac- 
ity variation with temperature, density and composition 
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with t he metal abundance usu ally scaled to solar compo- 
sitions. Ilglesias fe Rogers! lll993l) took a step forward for the 
OPAL project team by providing tables that include mix- 
tures enhanced in carbon and oxygen (C and O) relative to 
the base solar composition. We h ave incorporated their full 
range of composition tables in the lEggletonl JlQTlfj evolution 
code. In the last implementation of this code iPols et alJ 
Il995l) there only 10 of the 265 tables available were used. 
This provides us with the opportunity to refine our models 
so as to accurately follow the changes of opacity which occur 
in the later stages of evolution. The effects are expected in 
many types of stars. While these may be small for the evo- 
lution of main-sequence stars and white dwarfs we expect 
larger differences to be found in AGB and Wolf-Rayet stars. 

• Asymptotic Giant Branch (AGB) stars undergo third 
dredge-up which mixes helium burning products to the sur- 
face and forms carbon stars. Using the enhanced mixture 
tables we shall be able to model the thermal pulses and en- 
velope evolution more accurately. 

• Wolf-Rayet (WR) stars are massive and have lost their 
hydrogen envelopes exposing the helium cores. As time pro- 
gresses helium burning products are slowly exposed at the 
surface and in some cases the stars are eventually mostly 
composed of carbon and oxygen. 

We present our method of opacity table construction 
and detail its implementation in the Eggleton stellar evo- 
lution code. We discuss the effects on main-sequence stars, 
red giants and white dwarfs. We then present three tests of 
our work. The first is the effect of including extra carbon 
and oxygen on the structure and evolution of a low-mass 
population-Ill star, the second is a 5Mq thermally pulsing 
AGB star and the third a Wolf-Rayet star of 40M Q with 
mass loss. 



2 TABLE CONSTRUCTION. 

We start with the OPAL tables Jlglesias fc Rogerdll996l) as 
the framework around which we construc t our full tables i n 
an approach similar to that employed bv lPols et alJ lll995l) . 
However we choose to use the variable 72. = p/Tg rather than 
the density. This has the advantage that our tables can be 
a third smaller in this dimension. This partly compensates 
for increased memory requirement of the full tables. 

The OPAL tables extend in log 10 (T/K) from 3.75 to 
8.70 in steps of 0.05 to 0.2 and in log 10 (72/gcm~ 3 K~ 3 ) from 
— 7 to +1 in steps of 0.5. We choose the range of the final 
tables to cover from 3 to 10 in log 10 (r/K) with increments 
of 0.05, and from -8 to +7 in log 10 (72/g cm~ 3 K~ 3 ) with 
increments of 0.5. This gives us tables of 141 by 31 elements. 

To account for the effect of composition, the OPAL ta- 
bles include typically 265 T and 1Z grids for each metallicity. 
These tables are split into 5 groups with different hydrogen 
mass fractions X of 0.0, 0.03, 0.1, 0.35 and 0.7. The tables 
in each group have different compositions in helium, carbon 
and oxygen. 

The lowest temperature in the OPAL tables 
is log 1f i(r/K) = 3.75 a nd we use the tables of 
lAlexander fc Ferguson! (^394), which extend from 3 to 4.1 in 
log 10 (T/K), to complete them at low te mperature. Where 
the tables overlap they match well ilglesias fc Rogers! 



1996). It should be noted t hat the tables provided by 
lAlexander fc Ferguson! ill 994!) do not include enhanced 
carbon and oxygen mixtures. Therefore at temperatures 
when log 10 (T/K) < 4.0 our tables do not follow these 
mixtures. This is only important in the surface of AGB 
stars. We plan in future to deal with lo w temperature 
enhanced mixtures by including the work of iMarigq (l2002f) 
before application to the structure of AGB star envelopes. 

To complete each table we fill the region from 8.70 to 
10 in log i r| (r/K), where electron scattering dominates, ac- 
cording to lBuchler fc Yuen! l|l976f) . We then combine with a 
full table of the effective opacit y owing to electron conduc- 
tion for which we use the fits ofllbenl 111 975) to the tables of 
iHubbard fc Lamp! dl969T) and lCanutd \l97(t) . We combine 
the radiative opacity with the effective conductive opacity 
by reciprocal addition, 



^cff ^rad ^cond 

Before the complete tables are implemented in the stel- 
lar evolution code we perform one last process. We rescale 
the tables over the C and O plane to a more regular grid. 
In the OPAL grid the spacing used requires interpolation 
within irregular polygons, so to make the interpolation sim- 
pler and faster, we include extra tables in a similar but fully 
rectangular grid system. This means we have 61 tables for 
each X abundance or 305 tables for each metallicity. 



3 IMPLEMENTATION INTO A STELLAR 
EVOLUTION CODE. 

We impleme nt the tabl e s in t he stellar evolution code first 
de scribed bvlE g gletonl Jl97lfl and most recently updated 
by iPols et alTl]^95ll . Further details can be found in the 
references therein. The main drawback with implementa- 
tion of the new tables is the increased memory requirement. 
Whereas before there were only 10 opacity tables there are 
now 305. The new tables are each a third of the size of the 
old tables so the overall increase is a factor of ten. This 
corresponds to a memory requirement increase from 17 Mb 
to 170 Mb but the average specification of today's desktop 
computers comfortably accommodates this. One of the rea- 
sons why this memory requirement is so large is that we also 
store the spline coefficients for each of the tables. When we 
use the tables within the code we interpolate in the 1Z and 
T plane with bicubic splines. This ensures the smoothness 
of our tables. We then use linear interpolation in the 3 com- 
position dimensions to derive the final opacity value at each 
point. 

To test the effect of these tables we use three different 
interpolation schemes, 

• Method A uses 5 tables with X = 0, 0.03, 0.1, 0.35 and 
0.7 and Y — 1 — X — Z. Where X is the hydrogen mass 
fraction, Y the helium mass fraction and Z the initial, solar 
mixture, metallicity. We add a further two tables, both with 
X = Y = 0, the first with X c = 1 - Z and X = 0. The 
other with Xc = and Xo = 1 — Z. Where Xc and Xo are 
the mass fractions for the enhanced amount of carbon and 
oxygen above that included in the metallicity mass fraction 
so the total carbon mass fraction is Xc + Z(C) and for oxy- 
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gen Xo + Z(0) . This gives a total of 7 tables allowing simple 
interpolation in composition similar to the old method. 

• Method B is similar to method A but we include two 
extra tables for each hydrogen mass fraction. Again both 
have Y = 0. One has Ac — 1 — X — Z and Xo = and the 
other Ac = and Ao = 1 — A — Z. This leads to a total of 
15 tables in all. 

• Method C incorporates all 305 tables in the opacity cal- 
culations. We have 61 tables for each hydrogen abundance. 
When interpolating in composition we use three variables 
X, A c /(1 — X — Z) and A /(l - X — Z). The hydro- 
gen abundances are as for method A. We take values of 
0, 0.01, 0.03, 0.1, 0.2, 0.4, 0.6 and 1 for the carbon and oxy- 
gen grid planes. On this grid we perform a three-dimensional 
linear interpolation. 



4 TESTING & RESULTS. 

The new tables give almost identical results to the old for 
main-sequence stars and red giants. The timescale of core 
helium burning is altered. From method A to method C it 
increases by a factor of 0.1 percent for a 5M0 star. This 
is reassuring. Current stellar evolution models describe well 
the evolution of the main sequence and the red giant branch. 
However we do find that the opacity tables enhance numeri- 
cal stability in the later stages of evolution because they re- 
solve changes with composition in greater detail. This leads 
to a smoother variation in opacity. 

A similar result is found for white dwarfs. Between 
methods A and C for a carbon and oxygen white dwarf the 
radius increases by 0.2 percent. While conditions in the at- 
mosphere do pass through regions of the greatest difference 
between the two methods (the dotted line in figure [3J| there 
is not enough mass at these temperatures and densities to 
make a significant global difference. 

4.1 Polluting a O.5M Population-Ill Star. 

For our first detailed test we take a zero metallicity 0.5 Mq 
star and pollute it uniformly with carbon and oxygen. Such 
a star could have been formed just after the death of the 
first stars if the products of these stars were mainly carbon 
and oxygen. We choose such a low-mass star to limit the 
effect on the results of burning via the CNO cycle. Our ini- 
tial model has the following properties, A = 0.7, Y — 0.3, 
R = 0.45 R , T cff = 4780 K and L = 0.96 L . The tables 
provided by OPAL do not include one with A = 0.75, we 
could extrapolate to obtain the opacity but for this test we 
prefer to remain within the tables. From the initial model 
we perform three sub-tests. In the first we pollute the star 
at a constant C/O ratio of 1, then evolve each model and 
compare the time from zero-age main sequence to the helium 
flash for methods A, B and C. We also compare a sequence 
of models for stars of the same mass but with a metallicity 
of the same mass fraction as the carbon and oxygen. These 
results are presented in table Q We see that with the pol- 
lution fraction increase, the timescale of hydrogen burning 
increases because the energy produced at the centre takes 
longer to reach the surface. However when comparing the 
difference between methods B & C we see only a compar- 
atively small increase in timescale. When we compare with 



the scaled solar models we see a greater difference. This is 
due to the much larger opacity of heavier elements such as 
iron. 

With the second set of models we look at the radius, 
surface temperature and luminosity of the star on the main 
sequence as we pollute it with either only carbon or only oxy- 
gen (tables|2]and|2Jl. We use method A as the base to which 
we compare methods B and C. A similar trend is seen in the 
results in tables |2] and |H] We therefore see it is important 
to include the affect of enhanced C and O mixtures before 
hydrogen exhaustion. Similar mixtures are encountered in 
the envelopes of carbon stars of very low metallicity. 

Finally we consider the effect of varying the C/O ratio of 
the pollution (table [IJ at a constant pollutant mass fraction 
of 0.001. Table 2] demonstrates the need to use method C 
to follow accurately the variation in opacity as the C/O 
ratio changes. This is important in AGB evolution during 
the thermal pulses when the C/O ratio changes gradually as 
carbon becomes dominant. The difference between methods 
B and C in this case comes about because of non-linear 
structure in the C and O opacities that can only be resolved 
with the full set of tables. 



4.2 Thermal Pulses of an AGB Star. 

Our thermally pulsing AGB model is of a 5M star. To fol- 
low a number of pulses we includ e convective overshooting 
jSchroder. Pols fc Eggletonlll997l) during the pre- AGB evo- 
lution. This leaves the star with a larger core which reduces 
the rate at which the thermal pulses grow in strength. This 
means we do not find dredge-up in our models but we can 
evolve through the pulses. We obtain the average pulse pe- 
riod as the arithmetic mean of the times between the first 
ten helium shell flashes. Results are recorded in table 

First we find that the different methods affect the age of 
the star when it starts to undergo thermal pulses. Second we 
find that methods A and B have similar interpulse periods 
while method C's is 3.6 percent less. We attribute this to the 
change in opacity due to varying He/C/O ratios in the inter- 
shell region which alters the time evolution of the pulse. 

4.3 Evolution of a Wolf-Rayet Star. 

We have used the mass-loss rates of lNugis fc Lamersl (l2000t) 
to manufacture a Wolf-Rayet star of solar metallicity and 
an initial mass of 4OM . Some details are described by 
iDrav fc Tout! J2003). Our tables are expected to make a dif- 
ference because Wolf-Rayet stars become naked helium stars 
and helium burning products are mixed throughout and ex- 
posed at the surface of the star. We only use methods A and 
C and not method B because the carbon and oxygen mass 
fractions do not rise above those included in the base metal- 
licity until after the hydrogen envelope has been removed 
and so our tables only begin to make a difference after this. 

In figure we present the variation of opacity through 
the star, just after core carbon ignition, with mesh point, the 
grid over which we solve the equations of stellar structure. 
The solid line is with method A. We then take this struc- 
ture and use method C to calculate the opacities to give the 
dashed line. This model is not in hydrostatic equilibrium 
and the dot-dashed line shows the relaxed new structure 
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in hydrostatic equilibrium from method C. Comparing the 
physical parameters of these calculations we find the main 
change to be an increase in radius of 3.0 percent. The surface 
temperature has decreased by 1.5 percent. 

We then made two complete runs with similar initial 
conditions by methods A and C. In figure [5] we display the 
difference between Wolf-Rayet evolution as the hydrogen en- 
velope is removed. The greatest changes are in the radius 
and surface temperature. The star is in the WN phase be- 
tween 14 to 12Mq , after which it enters the WC phase. In the 
WN phase there are only small differences but upon entering 
the WC phase the differences build up so that the temper- 
ature is lower by around 2 percent and the radius about 1 
percent larger. The luminosity also tends to be larger by 
around 0.5 percent. Figure [3] is a plot of where the dif- 
ference between methods A and C are greatest in p and T 
when X = 0. The plot is of the mean difference in the opac- 
ity calculated by methods A and C at all possible carbon and 
oxygen abundances divided by the opacity when Y = 1 — Z. 
From the plot we see that the regions of greatest difference 
are small. The lines show where the Wolf-Rayet model cal- 
culated by method C lies at certain points of its evolution. 
The solid line is for the star at 12. 3M© and the dashes at 
6.IM0. From figure we see that the largest changes begin 
around the 12M0 point. This is when the model has moved 
into the region where the differences between A and C are 
greatest. 

We find the opacity tables alter the timescales for the 
late burning stages. Method C increases the helium burning 
timescale by 0.1 percent while the carbon-burning time scale 
increases by 2.5 percent relative to method A. These changes 
come about because the opacity varies with composition and 
effects the temperature structure and burning rates. While 
there is no large difference in a single variable the overall 
effect on structure and evolution is more significant. 
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5 CONCLUSIONS. 

Our major conclusion is that the changes induced by prop- 
erly including opacities for varying C and O mixtures are 
small. However they are a good thing to include in stellar 
evolution codes because they help numerical stability by re- 
moving sharp changes in the variation of opacity with com- 
position that are encountered when we use fewer tables and 
add another level of detail to models without much loss of 
computational speed. 

The main computational cost is the extra memory to 
store the tables. Method C requires ten times the memory 
(about 200 Mb) compared to methods A and B. However 
there is little computational speed cost: relative to method 
A, methods B and C require about 4 percent more time 
for the evolution of a Wolf-Rayet star even though we must 
evaluate four times as many opacities for each 1Z and T value 
than with method A. 

The tables and stellar evolution code are freely available 
from http:/ /www.ast.cam.ac.uk/^stars Interested users are 
welcome to contact the authors for details on how to down- 
load and implement our tables. 
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Figure 1. The variation of opacity throughout the Wolf-Rayet 
star just after core carbon ignition. The solid line is calculation 
from model A, the dashed line from model C and the dash-dotted 
lines from model C with hydrostatic equilibrium. The dotted line 
is the mass interior to the mesh point. 
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Table 2. Variation in physical parameters relative to method A 
with increasing carbon mass fraction. 
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Table 3. Variation in physical parameters relative to method A 
with increasing oxygen mass fraction. 
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Table 4. Variation in physical parameters relative to method A 
with varying the C/O ratio. 





AR/R 


AT efI /T eff 


AL/L 




/10~ 2 


/10~ 2 


/10~ 2 


10~ 6 


0.000 


0.000 


0.000 


IO" 4 


0.009 


0.023 


0.113 


10~ 3 


0.108 


0.209 


1.054 


ur 2 


1.140 


1.692 


8.721 


10~ 6 


0.000 


0.000 


0.000 


IO" 4 


0.009 


0.021 


0.108 


10~ 3 


0.106 


0.200 


1.013 


icr 2 


1.179 


1.710 


8.859 


10~ 6 


0.0000 


0.0000 


0.0000 


IO" 4 


0.0000 


0.0023 


0.0046 


10~ 3 


0.0023 


0.0092 


0.0410 


ur 2 


-0.0387 


-0.0181 


-0.1386 



Figure 2. The percentage difference in model parameters be- 
tween methods A and C versus mass of the Wolf-Rayet star. 
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Table 1. Variation of Evolutionary Time Scale with increasing 
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Table 5. Variation of timescales for a 5Mq star. 

Age at 1st Average Pulse 
Method Pulse / 10 s yrs Period / yrs 

A 1.199679 2 803 

B 1.193637 2 802 

C 1.200724 2 702 
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Figure 3. The value plotted by the contours is J S ( K C ~ k a) dXc dXo/ kq dXc (LXq. Here ka and Krj are the opacities calculated 
by methods A and C. The opacities are all calculated when X = and over the surface (S) defined by {Xq + Xq) = to 1. The vertical 
artifacts are from the wide spacing of the tables in p. We show two Wolf-Rayet models on this plot, solid line 12Mq and dotted 6.1Mq, 
in the same stage of evolution as in figure We see that the largest differences in the models indeed occur as the star drifts into the 
region where the opacity difference is greatest. Also plotted is the dotted line for a 0.3Mq pure C/O white dwarf model. While the line 
passes through regions of the diagram with greater values the white dwarf is only slightly affected because very little mass is in the 
regions with greatest difference. 



